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Abstract 

We present an improved measurement of the standard model CP violation 
parameter sin20i (also known as sin 2/3) based on a sample of 45 x 10 6 BB 
pairs collected at the T(45*) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. One neutral B meson is reconstructed in 
a J/ipK%, if){2S)K° s , XdK° s , r] c K° s , J/i;K* , or J/^K^ CP-eigenstate decay 
(y-j | channel and the flavor of accompanying B meson is identified from its decay 

■ products. From the asymmetry in the distribution of the time intervals between 

the two B meson decay points, we obtain sin20i = 0.82±0.12(stat)±0.05(syst). 

o 

CN 1 In the Standard Model (SM), CP violation arises from an irreducible complex 

phase in the weak interaction quark-mixing matrix (CKM matrix) |lj]. In particular, 
the SM predicts a CP-violating asymmetry in the time-dependent rates for B° and 
B° decays to a common CP eigenstate, fcp, with negligible corrections from strong 
interactions [01: 
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where T(B°,B° — > f CP ) is the decay rate for a B° or B° to fcp dominated by a 
b — > ccs transition at a proper time t after production, is the CP eigenvalue 
of fcp, Amj is the mass difference between the two B° mass eigenstates, and 0i 
is one of the three interior angles of the CKM unitarity triangle, defined as <pi = 
n — arg(— V-UVtd/ — V* b V c d). Non-zero values for sin20i were reported by the Belle 
and BaBar groups]^, |]]. 

Belle's published measurement of sin 20! is based on a 29.1 fb _1 data sample 
containing 31.3 x 10 6 BB pairs produced at the T(AS) resonance. In this paper, we 
report an improved measurement that uses 45 x 10 6 BB pairs (42 fb _1 ). The data 
were collected with the Belle detector at the KEKB asymmetric collider ||, which 

*Talk presented at Flavor Physics and CP Violation (FPCP), Philadelphia, U.S.A., May 16-18, 
2002. 
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collides 8.0 GeV e~ on 3.5 GeV e + at a small (±11 mrad) crossing angle. We use 
events where one of the B mesons decays to fcp at time top, and the other decays 
to a self-tagging state, /tag, i-e., a final state that distinguishes B° and B°, at time 
t tag . The CP violation manifests itself as an asymmetry A (At), where At is the 
proper time interval between the two decays: At = t C p — ttag- At KEKB, the T(4S) 
resonance is produced with a boost of (3^ = 0.425 nearly along the electron beam 
direction (z direction), and At can be determined as At ~ Az/(^)c, where Az is 
the z distance between the fcp and /tag decay vertices, Az = zqp — z tag . The Az 
average value is approximately 200 /j,m. 

The Belle detector || is a large-solid-angle spectrometer that consists of a silicon 
vertex detector (SVD), a central drift chamber (CDC), an array of aerogel thresh- 
old Cerenkov counters (ACC), time-of-flight scintillation counters (TOF), and an 
electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a su- 
perconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return 
located outside of the coil is instrumented to detect K° L mesons and to identify muons 
(KLM). 

We reconstruct B° decays to the following CP eigenstates f\ J/ipK®, ip(2S)Kg, 
XciK° s , r] c K° s for £ f = -1 and J/tpK° L for £ f = +1. We also use B° -> J/ipK*° decays 
where K*° — > Kgir . Here the final state is a mixture of even and odd CP, depending 
on the relative orbital angular momentum of the J/ip and K*°. We find that the final 
state is primarily = +1; the £/ = -1 fraction is 0.19 ± 0.02(stat) ± 0.03(syst) [[?]]. 
For reconstructed B — > fcp candidates other than J/ipK®, we identify B decays 
using the energy difference AE = — E^ m and the beam-energy constrained 
mass Mbc = \J (^beam) 2 — ip^) 2 ^ where E^ m is the beam energy in the center-of- 
mass system (cms), and Eg * and p c ™ are the cms energy and momentum of the 
reconstructed B candidate, respectively. Figure p] (left) shows the Mb c distributions 
for all B° candidates except for B° — > J/ipK^ that have AE values in the signal 
region. Table [I] lists the numbers of observed candidates (iV rec ) . 

Candidate B° — > J/ipK® decays are selected by requiring ECL and/or KLM 
hit patterns that are consistent with the presence of a shower induced by a neutral 
hadron. The centroid of the shower is required to be in a 45° cone centered on the 
K\ direction that is inferred from two-body decay kinematics and the measured four- 
momentum of the J/ip. Figure [I] (right) shows the p c ™ distribution, calculated with 
the B° — > Jji\}K\ two-body decay hypothesis. The histograms are the results of a 
fit to the signal and background distributions. There are 767 entries in total in the 
0.20 < 

pcms ^ 0.45 GeV/ c signal region |; the fit indicates a signal purity of 60%. The 
reconstruction and selection criteria for all of fcp channels used in the measurement 

T Throughout this paper, when a decay mode is quoted, the inclusion of the charge conjugation 
mode is implied. 

■^When the K\ is identified with the ECL only, the signal region is defined to be 0.20 < p ™ < 
0.40 GeV/c. 



2 




Figure 1: The beam-energy constrained mass distribution for all decay modes com- 
bined other than J/ipK® (left). The p c ™ distribution for B° — > J/ipK^l candidates 
with the results of the fit (right). 

are described in more detail elsewhere ||. 

Leptons, charged pions, kaons, and A baryons that are not associated with a 
reconstructed CP eigenstate decay are used to identify the 6-flavor of the accompa- 
nying B meson: high momentum leptons from b —>■ c£~T7; lower momentum leptons 
from c — > s£ + u; charged kaons and A baryons from b —>■ c —>■ s; fast pions from 
B° — > D ( -*)~(n + , p + , af, etc.); and slow pions from D*~ — > D°ir~. Based on the mea- 
sured properties of these tracks, two parameters, q and r, are assigned to an event. 
The first, q, has the discrete values q = ±1 that is +1 (—1) when B ta . g is likely to be 
a B° (B°), and the parameter r is an event-by-event Monte-Carlo-determined flavor- 
tagging dilution factor that ranges from r = for no flavor discrimination to r = 1 for 
an unambiguous flavor assignment. It is used only to sort data into six intervals of r, 
according to flavor purity; the wrong-tag probabilities, wi (I = 1, 6), that are used in 
the final fit are determined directly from data. Samples of B° decays to exclusively 
reconstructed self-tagged channels are utilized to obtain wi using time-dependent B°- 
B° mixing oscillation: (N OF — N SF ) / (N OF + N SF ) = (l — 2wi) cos(Am rf At), where A^qf 
and Asf are the numbers of opposite and same flavor events. The total effective tag- 
ging efficiency is determined to be J2i=i fi(l — ^ w i) 2 = 0.270 ± 0.008(stat)^gQQg(syst), 
where is the event fraction for each r interval. 

The vertex position for the fcp decay is reconstructed using leptons from J/ip 
decays or kaons and pions from r\ c and that for / tag is obtained with well reconstructed 
tracks that are not assigned to fcp- Tracks that are consistent with coming from a 
Kg — > 7r + 7r~ decay are not used. Each vertex position is required to be consistent with 
a run-by-run-determined interaction region profile that is smeared in the r-0 plane 
by the B meson decay length. With these requirements, we are able to determine 
a vertex even with a single track; the fraction of single-track vertices is about 10% 
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for zap and 30% for z taug . The proper-time interval resolution function, R sig (At), is 
formed by convolving four components: the detector resolutions for zqp and z tag , the 
shift in the z tag vertex position due to secondary tracks originating from charmed 
particle decays, and smearing due to the kinematic approximation used to convert 
Az to At. A small component of broad outliers in the Az distribution, caused by mis- 
reconstruction, is represented by a Gaussian function. We determine ten resolution 
parameters from the data from fits to the neutral and charged B meson lifetimes || 
and obtain an average At resolution of ~ 1.56 ps (rms). The width of the outlier 
component is determined to be (36l|) ps; the fractional areas are (6±1) x 10~ 4 and 
(3.1 ± 0.4) x 10~ 2 for the multiple- and single-track cases, respectively. 

After flavor tagging and vertexing, we find 766 events with q = +1 flavor tags and 
784 events with q = —1. Figure |2| shows the observed At distributions for the = 
+1 (solid points) and q£f = — 1 (open points) event samples. The asymmetry between 
the two distributions demonstrates the violation of CP symmetry. We determine 
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Figure 2: At distributions for the events with q£f = +1 (solid points) and q£f = — 1 
(open points). The results of the global fit with sin20i = 0.82 are shown as solid and 
dashed curves, respectively. 

sin 20i from an unbinned maximum-likelihood fit to the observed At distributions. 
The probability density function (pdf) expected for the signal distribution is given 
by 

e -\At\/r B o 

P sig (At, q, w h i f ) = — [1 - q£ f {\ - 2w,)sin 2fa sin(Am d At)], (2) 

At b o 

where we fix the B° lifetime (r^o) and mass difference at their world average values ||. 
Each pdf is convolved with the appropriate R sig (At) to determine the likelihood value 
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for each event as a function of sin20i: 



Pi = (1 - /oi) / [/si g ^sig(At', q, w h £ f )R sig (At - At') 

+ (1 - /si g )nk g (At')^bk g (At - At')]dAt' + /oiPoi(At), (3) 

where / s ig is the signal probability calculated as a function of p c ™ for J/ i\)K\ and 
of AE and M^ c for other modes. Pbkg(At) is the pdf for combinatorial background 
events, which is modeled as a sum of exponential and prompt components. It is 
convolved with a sum of two Gaussians, Pbkg, which is regarded as a resolution 
function for the background. To account for a small number of events that give large 
At in both the signal and background, we introduce the pdf, P a \, and the fractional 
area, / Q i, of the outlier component. The only free parameter in the final fit is sin20i, 
which is determined by maximizing the likelihood function L = Yl^Pi, where the 
product is over all events. The result of the fit is 

sin 20x = 0.82 ± 0.12(stat) ± 0.05(syst). 

The sources of the systematic error are listed in Table |2|. The systematic error is 
dominated by uncertainties in the vertex reconstruction. Other significant contri- 
butions come from uncertainties in the wrong tag fractions, the resolution function 
parameters and the J/ipK^ background fraction. 

A number of checks on the measurement are performed. Table [| lists the re- 
sults obtained by applying the same analysis to various subsamples. All values are 
statistically consistent with each other. Figure |3](a), (b), and (c) show the raw 
asymmetries and the fit results for all modes combined, (cc)Kg, and J/ipK®, re- 
spectively. A fit to the non-CP eigenstate self-tagged modes B° — > D^~n + , D*~ p + 
and J/ipK*°(K + 7r~), where no asymmetry is expected, yields 0.05 ± 0.04(stat). Fig- 
ure |3](d) shows the raw asymmetry for these non-CP control samples. 

Finally we comment on the possibility of direct CP violation. The signal pdf for a 
neutral B meson decaying into a CP eigenstate (Eq. (0)) can be expressed in a more 
general form as 



P sig (At, q, w u if) = e ' At ' TB ° {i + g (i _ 2wi) p^V^CP sin{AmdAt) 



At b o 



2|A|(-£/)a CP 



lAP + l 
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2 -l 


+■ 


A 


2 + l 



-cos(Am d At)]}, (4) 



where A is a complex parameter that depends on both B°-B° mixing and on the 
amplitudes for B° and B° decay to a CP eigenstate. The parameter acp in the 
coefficient of sin(Am,iAt) is given by acp = — £//raA/|A| and is equal to sin20i in 
the SM. The presence of the cosine term (|A| 7^ 1) would indicate direct CP violation; 
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Figure 3: (a) The raw asymmetry for all modes combined. The asymmetry for J/ipK^ 
and J/ijjK*° is inverted to account for the opposite CP eigenvalue. The corresponding 
plots for (b) (cc)Kg, (c) J/ipK^, and (d) non-CP control samples are also shown. 
The curves are the results of the unbinned maximum likelihood fit applied separately 
to the individual data samples. 



6 



the value for sin20i reported above is determined with the assumption |A| = 1, as 
expected in the SM. In order to test this assumption, we also performed a fit using 
the above expression with acp and |A| as free parameters, keeping everything else the 
same. We obtain |A| = LOlt^stat) and a C p = 0.82 ± 0.12 (stat) for all CP modes 
combined. This result confirms the assumption used in our analysis. 
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Mode 




□Kg 


J/^(£ + r)i^(7r+7r-) 

J/iP(£ + l-)K° s (ir°ir°) 

^(2S)(£ + £~)K° s (n + n-) 

jJj(2S) (J/^lT + lT-)K s (lT + lT-) 

ri c (K + K-ir°)K° s (n + iT-) 
rj c (K s K-ir + )K° s {n + iT-) 


636 
102 
49 
57 
34 
39 
33 


31.2 
20.8 
2.4 
4.3 
2.3 
11.1 
8.9 


J/^{£ + £-)K*°{K° s 7i°) 


55 


6.0 


J/^{£ + £-)Kt 


767 


307 



Table 1: The numbers of observed candidates (N rcc ) and the estimated background 
(iV bk g) in the signal region for each f C p mode. 



source 


+error 


—error 


vertex reconstruction 


+0.030 


-0.030 


wrong tag fraction 


+0.024 


-0.026 


resolution function 


+0.022 


-0.019 


background fraction (J/vfjK®) 


+0.014 


-0.015 


background fraction (except for J/ipK®) 


+0.007 


-0.006 


t b o and Arrid 


+0.007 


-0.006 


total 


+0.048 


-0.048 



Table 2: List of systematic errors on sin20i. 



Sample 


sin 20! 


= 5 U (q = +1) 
/ tag = B° (q = -1) 

J/lpK° s (7l + 7l-) 

(cc)K° s except J/ipK%(n + n-) 
J/^K° L 

J/i;K* (K s 7r ) 


0.60 + 0.19 
0.99 + 0.16 
0.67 + 0.18 
0.88 + 0.31 
1.14 + 0.23 
1.62 + 1.10 


All 


0.82 + 0.12 



Table 3: The values of sin20i for various subsamples (statistical errors only). 
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